LAMMPS (18 Feb 2013)
# Pour 2d granular particles into container

dimension	2
atom_style	sphere
boundary	f fm p
newton		off
communicate	single vel yes

region		reg block 0 100 0 50 -0.5 0.5 units box
create_box	1 reg
Created orthogonal box = (0 0 -0.5) to (100 50 0.5)
  1 by 1 by 1 MPI processor grid

neighbor	0.2 bin
neigh_modify	delay 0

# IMPORTANT NOTE: these values of stiffness (4000) and timestep (0.001)
# are used in this example file to produce a quick simulation and movie.
# More appropriate values for realistic simulations would be
# k = 4.0e5 and dt = 0.0001, as in bench/in.chute (if it were Hertzian).

pair_style      gran/hertz/history 4000.0 NULL 100.0 NULL 0.5 0
pair_coeff	* *

timestep	0.001

fix             1 all nve/sphere
fix		2 all gravity 1.0 spherical 0.0 -180.0
fix		xwalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 		xplane 0 100
fix		ywalls all wall/gran 4000.0 NULL 100.0 NULL 0 1 		yplane 0 NULL

region          slab block 1.0 99.0 30 34.5 -0.5 0.5 units box
fix             ins all pour 1000 1 4767548 vol 0.4 10 	        diam 0.5 1.0 region slab
Particle insertion: 224 every 3000 steps, 1000 by step 12001

fix             3 all enforce2d

compute		1 all erotate/sphere
thermo_style	custom step atoms ke c_1 vol
thermo		1000
thermo_modify	lost ignore norm no
compute_modify	thermo_temp dynamic yes

#dump		id all atom 250 dump.pour

#dump		1 all image 500 image.*.jpg type type #		zoom 1.6 adiam 1.5
#dump_modify	1 pad 5

run		25000
Memory usage per processor = 9.58 Mbytes
Step Atoms KinEng 1 Volume 
       0        0           -0            0         5000 
    1000      224    185.87695            0         5000 
    2000      224    348.08417            0         5000 
    3000      224     565.2162            0         5000 
    4000      448    1042.9716            0         5000 
    5000      448    1542.0614            0         5000 
    6000      448    1766.5136    41.192729         5000 
    7000      672    1846.4596    67.352592         5000 
    8000      672    1814.3429    92.898167         5000 
    9000      672    1662.3643    111.88767         5000 
   10000      896    1654.7184    91.390546         5000 
   11000      896    1532.7703    69.886692         5000 
   12000      896    1470.9688    66.246593         5000 
   13000     1000    1471.1825    60.129537         5000 
   14000     1000    1355.6479    59.801716         5000 
   15000     1000    1132.1071    52.624027         5000 
   16000     1000    773.40416    42.626873         5000 
   17000     1000    461.51805    37.889877         5000 
   18000     1000     338.7731    28.037959         5000 
   19000     1000    166.73672    26.940654         5000 
   20000     1000    45.061374    18.328659         5000 
   21000     1000    24.350861    10.312694         5000 
   22000     1000    18.811191    8.8044879         5000 
   23000     1000    15.483854    5.3046832         5000 
   24000     1000    11.962946    5.0549062         5000 
   25000     1000    9.0369801    3.2724609         5000 
Loop time of 2.88655 on 1 procs for 25000 steps with 1000 atoms

Pair  time (%) = 1.22176 (42.3259)
Neigh time (%) = 0.303655 (10.5196)
Comm  time (%) = 0.0110941 (0.384337)
Outpt time (%) = 0.000306368 (0.0106136)
Other time (%) = 1.34974 (46.7595)

Nlocal:    1000 ave 1000 max 1000 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Nghost:    0 ave 0 max 0 min
Histogram: 1 0 0 0 0 0 0 0 0 0
Neighs:    2330 ave 2330 max 2330 min
Histogram: 1 0 0 0 0 0 0 0 0 0

Total # of neighbors = 2330
Ave neighs/atom = 2.33
Neighbor list builds = 2001
Dangerous builds = 0
